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ABSTRACT 

Future high spectral resolution telescopes will enable us to place tight constraints on 
turbulence in the intra-cluster medium through the line widths of strong emission 
lines. At the same time, these bright lines are the most prone to be optically thick. 
This requires us to separate the effects of resonant scattering from turbulence, both 
of which could broaden the lines. How this can be achieved has yet not been quanti- 
tatively addressed. In this paper, we propose a flexible new parametrization for the 
line profile, which allows these effects to be distinguished. The model has only 3 free 
parameters, which we calibrate with Monte-Carlo radiative transfer simulations. We 
provide fitting functions and tables that allow the results of these calculations to be 
easily incorporated into a fast spectral fitting package. In a mock spectral fit, we ex- 
plicitly show that this parameterization allows us to correctly estimate the turbulent 
amplitude and metallicity of a cluster such as Perseus, which would otherwise give 
significantly biased results. We also show how the physical origin of the line shape can 
be understood analytically. 
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1 INTRODUCTION 

Turbulence plays a key role in many aspects of the 
physics of the intra-cluster medium (ICM). It pro- 
vides non-thermal pressur e support, thus biasin g the 



Despite its importance, turbulence in the ICM has 
so far only been weakly constrained, generally by in- 
direct means. These have came from the analysis of 



masses measured in X-ray (|Lau. Kravtsov fc Naga. 
and via the S unyaev -Z eldovich effect dShaw et al 
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iBattaglia et~ai1 l2011al fbl iParrish et ail |2012| ). which 
sume thermal pressure alone. This uncertainty in cluster 
mass calibration constrains the use of clusters for pre- 
cision cosmology. Turbulence could also transport met- 



the lack of detection of res onant- line scattering 

20lbT). Faraday ro - 



( Churazov et al. | 12004: IW crncr ct al 



tation maps i|Vogt fc Enfflinl 120051 : lEnfilin fc Vogtl l2006h . 



als l|Rebusco et al.l 120051; ISimionescu et al.l l200Sf). amplify 
magnetic fields dSubramanian.' Shukurov fc Haugenl [20061; 



iRvu et al.1 [20081 : ICho et all [20091: iRuszkowski et al.ll201ll ). 

and by accelerating cosmic rays, is likely to be closely 
related to radio phenomena i n clusters (|Brunetti et al.l 
l200ll : iBrunetti fc Lazarian|[2007l ) . Moreover, turbulence can 
play a crucial role in the thermal state of the cluster. 
It can h alt catastrophic cooli ng by dissip ation of tur- 



bulence (|Churazov et al. 20041: llCunz et al.l I2011D. turbu- 

' Kim fc Naravan|[2003: 

more subtly, by randomiz- 



lent heat diffusio n, rtCho et al.H2003l 
m 



iDennis fc Chandranl 12005^ or 
ing the magnetic field, and resto ring thermal conduction 
to ~ 1/3 of the Spitzer value (Ruszkowski fc Ob] l20ld. 



and deviations from hydrostatic equilibrium with therma l 
pressure alone ([Churazov et al. 1 120081 : IChurazov et~al]|2010l : 
IZhang et ai1l200Sh . Constraints can also be provided by a 
more direct means - measuring the Doppler broadening of 
strong emission lines. Unfortunately, due to the poor energy 
resolution of current X-ray telescopes, to date this type of 
analysis has only yielded upper limits (|Sanders et alj|2010l : 
ISanders. Fabian fc Smith! l201ll : Ide Plaa et al.l |2012| ). The 
challenge stems also from the fact that the expected gas mo- 
tions are highly subsonic, with Mach numbers A4 ~ 0.1 — 0.3 
in the cluster core. 

However, this situation is expected to change with the 
launch of Astro-H3 in 2014 and ATHENA3 in the farther 
future. With their high spectral resolution (e.g., FWHM ~ 



120111 : IParrish. Quataert fc Sharmall2010l ). Amongst others, 
gas motions can be seeded by structure formation (including 
motions of sub halos and mergers) and AGN activity. 
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4 — 5 eV for Astro-H), these telescopes can accurately mea- 
sure the widths of emission lines, thereby placing tight con- 
straints on the amplitude of turbulence. They may even al- 
low us to separate different modes of g as motions by mixtur e 
modeling of the detailed line profile |Shang fc Oh|[2012al V 
For instance, bulk motions and gas sloshing can be identified 
through components with different means, while the volume 
filling fraction of turbulence can be identified through com- 
ponents of different widths. To some extent, the high spec- 
tral resolution of Astro-H can mitigate against its poor spa- 
tial resolution in learning about the detailed velocity field. 

To correctly recover the amplitude of turbulence from 
emission lines, one important factor has to be taken into 
account: resonant scattering. The ICM is optically to thin 
continuum free-free emission. However, the optical depths 
for a few strong metal emission lines could still exceed unity 
urazov et al.ll201(il . and references therein). Photons in 
the lines therefore experience multiple scatterings before 
they reach the edge of the cluster. Such RS has many impor- 
tant consequences. For instance, the line ratios of optically 
thick and thin lines are altered. The line width and shape are 
changed as well, since photons in the optically thick line cen- 
ter are scattered to the optically thin line wings, where they 
can escape. Morever, these effects happen systematicaly as a 
function of position, since photons escape more easily from 
the cluster outskirts than from the center. Under certain 
conditions, the variance of the line could increase by 40 %. 
Failure to account for resonance scatterin g introduces biases 
in th e measurement of the metallicity l|Sanders fe Fabian] 
2006) and turbulence (as we shall show). 

RS in the context of galaxy clusters has been 
studied both analytically and with radiative-transfer 
Monte-Carlo simulations. It was fi rst considered by 
iGilfanov. Svunvaev fc Churazovl 0987), w ho discussed var- 
ious featur es of the RS effect and their astrop hysical im- 
plications. Mathews . Buote fc Brighentil (2001) considered 
if the RS could suffi ciently explain the abundance hole o f 
metals seen in M87. ISazonov. Churazov fc Sunvaevl (|2002T ) 
demonstrated that Rayleigh scattering could polarize the 
line emission, and detection of the polarized X-rays could 
yield valuable information on ICM physical conditions . This 
work was later followed up bv lZhuravleva et alj ((2010) , who 
proposed to use the polarization of X-ray lines as a way of 
measuring the tangential compon e nt of the ICM velocity. 
iMolnar. Birkinshaw fc Mushotzkvl (|2006T ) discussed how the 
resonant scattering lines could be used to measure angular 
diameter distances to clusters and consequently constrain 
cosmology. ISanders fc Fabian] (I2006T ) considered the effect of 
RS on the inferred metallicity of the cluster, via an analytic 
calculation which can used to fit observations. This work is 
most similar in spirit to this paper, although they adopted 
the single scattering approximation, ignored gas motions, 
and did not include the l i ne bro adening due to RS. Most 
recently, IZhuravleva et all |2011ai ) investigated how the RS 
depends on the spatial scale and anisotropy of gas motions 
and how such dependence could be used to probe the char- 
acteristics of gas motions. 

Despite the increasing sophistication of these studies, 
perhaps the most important leading order effect-the joint 
effects of turbulence and RS in altering the line profile — has 
not been quantitatively addressed. Studies have generally 
either considered the effects of RS in a static ICM, or the 



effect of gas motions in an optically thin ICM. In reality, 
the gas is likely to always have subsonic motions, and the 
brightest lines with the highest oscillator strengths (e.g., the 
He-like line of Fe XXV) we will use to measure these mo- 
tions are by the same token the most likely to be optically 
thick, and the effects of RS cannot be ignored. It is im- 
portant to consider these effects in tandem for two reasons. 
Firstly, as previously mentioned, by altering line widths, RS 
could bias measurements of turbulence. Our first purpose in 
this paper is to learn if one can strip away the confound- 
ing effects of RS to learn about the underlying velocity field. 
Secondly, the two effects interact with one another: turbu- 
lence reduces the line optical depth and hence the effects of 
RS, while RS alters the effective velocity shift experienced 
by photons. Our second purpose is therefore a practical one: 
we aim to provide observers a parametric model of the line 
profile which takes both RS and turbulence into account, in 
order to quickly and accurately estimate turbulent broaden- 
ing. The alternative would be to run radiative transfer sim- 
ulations for each plausible realization of a given cluster's 
velocity field, a very costly and slow way of fitting the data. 

The remainder of this paper is organized as follows. 
§ [2] describes the assumed cluster model, the line profile 
model and procedures of Monte-Carlo simulation and pro- 
file fitting. § [3] discusses how the line profile depends on 
various physical parameters. We then calibrate the depen- 
dence with a suite of simulations. The results are given in 
both tabulated and analytical fitting function form. We also 
test the model and the calibration by comparison with a 
new set of simulations, performing a mock spectral fitting. 
§ 13.61 describes how the physical origin of the characteristic 
line shape after resonant scattering can be understood ana- 
lytically. § [3] summarizes our results and discusses possible 
extensions to the work. 



2 METHODOLOGY 

In this section, we describe our method of a flexible paramet- 
ric model to describe the line profile in the presence of gas 
motions and resonant scattering. We assume that the clus- 
ter temperature, density and metallicitjjf] profiles are known 
ahead of time from X-ray observations, and only the velocity 
field is unknown. Given realizations of the velocity field mo- 
tivated by numerical simulations, we fit a line profile which 
separates the effects of turbulent broadening from that of 
RS. 

2.1 Cluster Model 

Electron density profile The isothermal f3 model has 
been widely used to describe th e free-free surface brightness 
and d ensity profiles of clusters (|Cavaliere fc Fusco-Femiand 

Eezi), 

/ "3/3/2 

n(r) = n (l + -J , (1) 

/ J\ 1/2-3/3 

S(r) = S (l + 7l j , (2) 

3 We justify this assumption later. 
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Here no and So are the central density and surface bright- 
ness, respectively, and r c is the core radius, while /3 con- 
trols the asymptotic density profile at large radii. The model 
is very simple with only three free parameters. Overall, it 
describes observed profiles reasonably well, although the 
density is often over(under)-predicted at the large(small) 
radii. For this reason, several extensions have been proposed, 
which fit obser ved profiles better at the cost of more fr ee pa- 
rameters (e.g. JChen et al.ll2007l ; IVikhlinin et~al1l2006l ). 

In this first attempt to systematically model the line 
profile with resonant scattering, we restrict ourselves to the 
simple P model. Since the line optical depth is only linearly 
proportional to density, we do not expect the results in this 
paper to be strongly affected by deviations from the j3 model. 

Ion density profile Resonant scattering (oc nzj and 
line emission (oc n^n^) is sensitive to the density of metal 
ions nz ; . This is sensitive to two factors: metallicity (nz oc 
Z), and temperature (which sets the ionic abundance in 
coronal equilibrium). In this paper, we shall also assume 
that nzj can be fit by a f3 profile (with different parame- 
ters from the electron profile). While not true in detail, it is 
reasonably accurate for the central core where most of the 
emission and resonant scattering takes place; the beta profile 
can be viewed simply as a mathematical fitting function for 
nzj with a central core and asymptotic outer slope, without 
reference to the electron /3 profile. In practice, nz ; oc n e is 
a good approximation for NCC clusters, as we shall verifjjf]. 
Since this is mathematically simpler, we treat this case first. 
We consider the NCC case where nzj , n c have different /3 
model parameters in M3.4I 

Temperature The temperature profile has two effects: 
affecting ionic abundance in coronal equilibrium (discussed 
above), and making thermal broadening a function of ra- 
dius. The former effect is the dominant one. For the latter, 
we initially assume spatially constant thermal broadening 
(i.e., isothermal temperature, which is approximately true 
in the core), and then show that taking the true tempera- 
ture profile into account has little effect. 

Velocity Field We generate realizations of turbulence 
over a 3D grid of 256 3 , with a box size of 2.2 Mpc. The 
velocity fields are generated from a power spectrum of the 
form: 



E(k) = 



No 



Mi)' 



k kc 



(3) 



This form is motivated by both simulations and indirect ob- 
servations of pressure fluctuations which fin d power-law like 



spectra and a cutoff in the large scale (e.g., Schuccke r et al 



120041 : IVazza et al1l201ll : IVazza. Roediger fc Brueggenll2012j r 

For Kolmogorov turbulence, the slope n s = —5/3. 

The use of a power spectrum alone to describe turbu- 
lence is likely to be incomplete and over-simplified, since it 
does not take possible spatial dependence of the gas mo- 
tions into account, particularly radial dependence as well 
as anisotropy (which is important since we can only probe 
gas motions along the line of sight). Strictly speaking, the 



4 This arises simply because the gas is approximately constant 
metallicity and isothermal in the core of NCC clusters, and the 
first order variations (declining metallicity and decreasing tem- 
perature with radius) tend to cancel out. 



velocity field in a cluster is likely to be non-Gaussian. A 
more realistic approach would be to take the velocity fields 
directly from numerical simulations. However, this simple 
prescription is in relatively good agreement with numerical 
simulations, and moreover allows us to rapidly parameter 
space. For relatively simple statistics such as line height and 
width, a power spectrum parametrization is adequate. This 
is less likely to be true for higher order statistic whic h de- 
pend on details of the line shape (|Shang fc Ohll2012bl ). 

Another possible concern lies in the assumption that 
turbulence does not exhibit systematic spatial trends in the 
cluster. This seems in contradiction with recent simulation 
results (e.g., lLau. Kravtsov fc Nagail [20091 ; iBattaglia et al.l 
l2011al ). which show a power law increase toward large radii. 
However, the power law increase usually takes place around 
0.5 — Irsoo (0.5-1 Mpc for nearby bright clusters), which is 
considerably larger than the core radius ( < 200 kpc). We 
therefore do not expect it to significantly alter our conclu- 
sions. 

The normalization factor No is adjusted to obtain the 
desired RMS velocity. From numerical simulations, we gen- 
erally expect turbulent gas motions to be highly subsonic, 
with Mach numbers M ~ 10 - 30%, although the Mach 
number rises as a f unction of radius and can be transonic 
at the virial radius ( Sha w et al.ll2010r ) . We also place an up- 
per bound on sub-grid turbulence by integrating the power 
spectrum from fc ma x (the wavenumber corresponding to the 
smallest cell; in our simulations the Nyquist length scale is 
8.6 kpc) to infinity. The velocity dispersion of this sub-grid 
turbulence, cr su b, will be used in the following steps. 



2.2 Line Profile Model 

Resonant scattering (RS) modifies the line profile by re- 
distributing photons from the optically thick line center to 
the optically thin line wings, where they can escape to the 
observer. Here, we describe our fitting function to the line 
shape. 

In the absence of resonant scattering, the emission line 
profile can be well approximated by a Gaussian, so long as 
the injection scale of turbulence is smaller than the emit- 
ting region where most photons arise (thus, the field of view 
encompasses many independent 'patches' and the central 
limit theorem holds). RS modifies the line shape due to its 
frequency-dependent cross section S (the symbol £ is used 
to distinguish it from the line dispersion a), 



Z(E,a) 



1 



2na 



(4) 



where Eo is the central energy of the line and a is the disper- 
sion of the line, including contributions from both thermal 
and turbulent motions: 



22 2 -^0 / kT 2 \ 

a = a ther + a turb = ^ ^turb) 



(5) 



Here, a" tur b is the ID RMS velocity dispersion, c is the 
speed of light, m p is the proton mass, and A is the atomic 
number. In the discussions that follow, we specialize to the 
iron ion (A = 56), in particular the He-like 6.7 keV line, 
which for instance in the absence of turbulent motions has 
a line center optical depth of 2.7 in Perseus and 1.4 in Virgo 
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|Churazov et alj|20ld ). However, our methods could be ap- 
plied to any other resonance line. Photons in the optically 
thick line center are scattered out of the line of sight, while 
photons in the optically thin line wings escape freely toward 
the observer. During the scattering process, photons which 
are re-emitted at line center can continue to scatter, while 
photons which are re-emitted in the lines escape. Thus, the 
line center is suppressed, while the wings are slightly en- 
hanced due to re-emission. This is clearly seen in Fig. \T\ 
We therefore model the line profile as: 



P(E, a) = A(E, a)R(E, a)G(E, a) 



(6) 



where G(E, a) is the original Gaussian line profile without 
RS, and A(E,a) and R(E,a) are factors accounting for ab- 
sorption and reemission, respectively. We found the modi- 
fication to the line profile due to re-emission could be ad- 
equately modeled with a simple power law with two free 
parameters: 

(7) 

where to is the opacity at line center from the cluster center 
to infinity, while £ = tq<t (since to oc E oc a -1 , we set 
ro = C/ ")- We separate a from tq in the equations, as it 
is the quantity of interest in this paper. £ can be readily 
computed in terms of fundamental constants and observed 
quantities: 



R(E, a) = (1 + 9l r h ) = (1 + 9i((/a) h ), 



a l £(£ = E )n,idl 



y/^hr e cfSi(l)n z (l)dl, (8) 



where h is the Planck constant, r e is the classical radius of 
the electron, / is the oscillator strength of the transition, 
Si is the fraction of the element in the appropriate ioniza- 
tion state, and nz is the number density of the element. 
We model the energy-dependent absorption using the usual 
exponential function with an additional free parameter, 



A(E,a) = exp(—g 2 T(E)) = exp — g 2 — e 2^ 

a 



(9) 



In total, there are three free parameters: g\,g 2 ,h. We find 
this model is sufficiently flexible to fit the simulated profiles. 

2.3 Monte Carlo Radiative Transfer 

For resonance line radiative transfer where multiple 
scatterings are possible, Monte Carlo simulations offer 
the greatest flexibility and ease of computatio ns, and 



Bonilha et al.l 



are standard i n the fi eld (e . . 
Ahn. Lee fc Led ll200ll ): I Zheng fe Miralda-Escudi 
Hansen fc ObTl|2006l )). T he simulations performe d for 



Till); 



2002); 



this work largely follow [Z huravlcv a et al 



also ISazonov. Churazov fc Sunvaevl ( |2002h : IChurazov et aLj 



l|2004l ): IZhuravleva et all (|2010l l). For completeness, we nev- 
ertheless describe briefly the main procedures. A simulation 
is composed of the following steps: 
1) Initialization of photons. 

We produce a mock photon with a unit weight w, pro- 
pogating in a random direction. Its location is randomly 
drawn according to the volume emissivity distribution, while 
its energy is taken from a Gaussian distribution with a vari- 



ance of (T th „ + cr t , 



The mean of the Gaussian distribution 



takes into account the Doppler shift due to cell motion, 

Emean = E ( 1 + V ce ll-m/c), (10) 



where v ce ii and m are the vectors of the cell velocity and 
photon propogation direction, respectively. We record the 
initial location and energy of the photon. This information 
is used to extract line profiles without RS, which can be 
compared to profiles with RS to determine model parame- 
ters. 

2) Photon propagation. 

Next, we integrate along the propagation direction m to 
find the optical depth to the edge of the cluster and the cor- 
responding escape probability p esc = e~ T (the cluster is cut 
off at 1 Mpc). The weight of the escaped photon is then equal 
to wp esc . Alternatively, the photon might be scattered. The 
optical depth to the next scattering location, T nex t, follows 
an exponential distribution between and r. We therefore 
randomly draw T next = — ln(l — £(1 — p eac )), where £ is a 
uniformly distributed random number between and 1. The 
position of the next scattering is then identified using T nex t 
and m. 

3) Scattering. 

RS takes place and changes the direction of photon 
motion and its energy. The new propogation direction of 
the scattered photon m' is randomly drawn according to 
the proper scattering phase matrix. There are two common 
types of scattering: isotropic and Rayleigh (dipole) scatter- 
ing. The scattered direction of the former is uniformly dis- 
tributed, while the probability of m' for the latter is pro- 
portional to (m • m') 2 . Their relative weights depend on the 
total angular momentum change j of the ground state and 
the difference Aj between the ground and excited state. For 
the He-like 6.7 keV like of iron, scattering is purely dipole. 
Thus, in this paper we mainly focus on Rayleigh scattering, 
while results for the isotropic scattering will be posted on 
the interne10. To find the energy of the scattered photon 
E' , we first need to compute the velocity of the scattering 
ion. Along the original propogation direction m, the veloc- 
ity of ion must satisfy Eo = E(l — v ion ^/c) in order for 
the RS to happen. Perpendicular to m, the ion velocity is 
drawn from a Gaussian distribution of mean v ce ii and vari- 
ance er 2 her + o"sub- The new energy is then computed using 
m' and Vi on , E' = Eo(l + v ce ii • m'). The scattered photon 
has a weight w' = w(l — p esc ). 

4) Repeat 

We repeat step 2-4 many times until the weight of the 
photon is neglible (< 10 -8 ). We run the simulation for a 
large number of photons (10 5 — 10 6 ), which allows us to 
build up good statistics to compute the PDFs and fit the 
model parameters. 



I 2011bl) (see 2.4 Fitting Procedure 



We fit the simulated line profiles to obtain the best-fit pa- 
rameters and their uncertainties. As these fits are performed 
many times, we choose a simple and fast strategy: we min- 
imize the x 2 between the model and binned data using the 
Levenberg-Marquardt algorithm implemented in the GSL 
scientific librarj|j. In particular, we minimize: 



= E 



(N air , 



N mod.i ) 



N a 



(11) 



5 https:/ /sites. google.com/site/cienshang/toolbox 
e http:/ /www. gnu.org/software/gsl/ 
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Figure 1. Line profiles with and without RS from the Monte 
Carlo simulation, together with our best fit model. The photons 
are accumulated from ring with a width of 10 kpc, centered at 
55 kpc. Please see text for other details; the analytic profile is 
discussed in § 13.61 



where N 3 im,i and N m od,i are simulated and model-predicted 
photon counts in i-th bin, respectively. Since x 2 minimiza- 
tion is inaccurate for small number st atistics (wh ere one 
should minimize the Cash C statistic l|Cashlll979t ). which 
correctly maximizes the Poisson likelihood), only bins with 
Nsim > 10 are used in the fit. 



G(E, a) unmodulated by the effects of resonant scattering 
can be recovered. 

In order to reduce the dimensionality of the problem, it 
is useful to exploit self-similarity. There are 4 length scales 
in the problem: the photon mean free path at line center 
from the cluster center A, the projected radius r, the field of 
view FOV, and the core size r c . The solution only depends 
on ratios of these length scales; for convenience, we choose 
to fix the core size r c , and vary A,FOV,r. For definiteness, 
we set r c — 100 kpc, though obviously our solutions are in- 
dependent of this choice, and would hold with appropriate 
rescaling. The photon mean free path A depends on the op- 
tical depth at line center To. Since to = C/c, we can choose 
to vary the overall normalization of the optical depth pro- 
file to by fixing £ and varying a. Finally, given a model of 
the radially varying line profile collected in an infinitesimal 
area P(r), the observed profile is given by integrating the 
emission- weighted profile P(r) over the FOV: 



P(FOV) = 



J POV S(r)P(r)d 2 r 
f FOV S(r)d*r ■ 



(12) 



Note P(r) is collected in a 2D ring on the sky, not in 3D 
shells. We check that self-similarity indeed holds numerically 
in |T3] 

The dependence on the remaining four parameters 
(n s ,k c ,/3 and r) are however not as simple, and need to be 
determined by Monte Carlo simulations. Note that /3 and r 
are known to us while n 3 and fc c are generally not known. We 
therefore discuss these two classes of parameters separately. 



3 RESULTS 

The main task of this paper is to present a new model for the 
line profile with RS, and to calibrate the dependence of the 
model on various physical quantities. We first consider the 
case when nz i oc n e , since this is mathematically simpler, 
and discuss the case where nz; has a separate profile later. 

The relevant quantities are: the cluster temperature 
profile T(r), f3 model parameters for the electron and ion 
profiles (no, j3, r c ) - note that the beta model parameters for 
the ions implicitly folds in information about the metallicity 
and temperature profiles), velocity field parameters (a tur t, 
n s , fc c ), and geometrical parameters parameters (radius r, 
and field-of-view (FOV) of the telescope). Dependence on 
some of the parameters has already been included in the 
model (equation [5] - O itself. For example, the dependence 
on temperature and aturb is approximately incorporated^ 
through a: both higher temperature and stronger turbulence 
raise a. Likewise, the normalization of the ion profile, m,o 
is set by £ (see equation ij8}. In the simulations, we fix T 
to 4 keV {(Jther ~ 1.8 eV for iron ions) and set £ = 5 eV 
(see discussion below) . We then fit the line profile for various 
values of aturb (and consequently a). If our proposed model 
accurately fits the line profile, and factors A(E , a) , R(E , a) 
appropriately describe the effects of resonant scattering and 
can be accurately recovered, then the original line profile 

7 This is exactly true if the velocity and temperature profiles are 
constant, and approximately true otherwise. See § 13.11 for more 
details. 



3.1 Dependence on n s and k c 

First, we run a series of simulations of 10 6 photons with dif- 
ferent k c . In these simulations, n s is fixed to its Kolmogorov 
value, ns = —5/3. We collect photons from a set of con- 
centric rings with width of 10 kpc. This choice of width 
is narrow compared to r c , and so can be approximated as 
"infinitesimal" (the goodness of this approximation will be 
checked in § 13 - 3ft , but still sufficiently wide to collect enough 
photons for good statistics. The inner edges of these rings 
are at r in = 0, 50, 100, 300 kpc. The result for r ln = 100 
kpc is shown in Fig. [2] which shows a transition around 
fc c ,crit ~ 2-K/r c = 0.0628 kpc" 1 . Before this transition, the 
model parameters vary with k c , while after it, they are con- 
sistent with constant. The same conclusion could be drawn 
from other cases. This result could be understood intuitively. 
At small k c , the outer stirring scale is larger than the size of 
the emitting region (~ r c ), and the LOS only "sees" part of 
an eddy. So the motion is more similar to ordered bulk flow 
than chaotic. The results in this regime depend on the de- 
tailed configuration of the velocity field, and each case has 
to be studied individually (in particular, a mixtu re model 
approach to fitting the PDF l|Shang fc Ohll2012bh may be 
applicable here). At larger k c , on the contrary, the eddy 
size is small relative to the size of emitting region; since 
photons transverse many eddies, by the central limit theo- 
rem, the assumption of random motions which have a ID 
gaussian profile is valid. In this regime, the line shape no 
longer depends on k c or o t her I 'aturb- Indeed, the model pa- 
rameters are also independent of n s . This is seen in Fig. [3] 
which shows the parameters as functions of n s , assuming 
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Figure 2. Model parameters as functions of k c . The photons are 
collected from a ring with a width of 10 kpc and inner edge at 
100 kpc. n s is fixed to -5/3. 
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Figure 4. Reduced \ 2 f°r au th e fits. 
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Figure 3. Model parameters as functions of n s . In this figure, k c 
is fixed to 0.1 kpc -1 . 



fc c = 0.1 kpc" 1 . Hereafter, we will focus only on this regime, 
and fix k c = 0.1 kpc -1 , and n s = —5/3. 

This relative independence to the details of the velocity 
power spectrum is easy to understand. Once k c > fc CiCr it, so 
that the velocity field sampled by outgoing photons is ap- 
proximately Gaussian, then the velocity field largely affects 
the PDF by changing its line width (note that all these cal- 
culations are performed at fixed line width, by adjusting the 
power spectrum normalization No; equation the PDF is 
not sensitive to the detailed spatial structure of the velocity 
field. If we were to consider higher-order statistics which are 
sensitive to the spatial st ructure of the velocity field, such 
as the structure function l|Zhuravleva et al|2012i ), we would 
not observe this behavior. In this case, the independence of 
the details of the velocity field is a great advantage - the ef- 
fects of resonant scattering would otherwise be much more 
difficult to compute. 



3.2 Dependence on f3 and r 

To determine the dependence on /3 and r, we run a suite of 
simulations, varying /3 from 0.5 to 1 with steps of 0.1, and r 
(the radii of the rings from where the photons are collected) 
from 5 kpc to 705 kpc with steps of 50 kpc. For a given set 
of p and r, we set the r.m.s. of the turbulent velocity to 0, 
100, 200, and 300 km/s, covering a range of to from 2.7 to 
0.7 (note that our model automatically recovers the original 
Gaussian profile as to — > 0, so the zero point does not need 
to be checked). Each set of /? and r has 4 simulated profiles, 
which we fit simultaneously to constrain pi, f>2, and h. The 
results are tabulated in Table [ATI and lA2l in the Appendix. 
An example is shown in Fig. [1] with f3 = 0.6, r = 55 kpc and 
o~turb = 0. Fig. [4] shows reduced \ 2 f° r au the fits. The fits 
are in general very goocjfl with the mean reduced \ 2 = 0.74 
Equation [6] - [9] 1121 and Table IA1I are everything one 
needs to build a line profile model for spectral fitting. For 
an arbitrary f5 and r/r c , interpolation needs to be performed 
to determine the proper pi, p2, and h. Alternatively, we pro- 
vide below polynomial fitting forms of pi(/3, r), P2(P, r) and 
h(f3,r), found by scientific data mining software Eureqa 



8 The somewhat better than expected reduced \ likely comes 
from the simulation scheme of assigning weight to photons, which 
may artificially suppress the shot noise in the simulated profile. 
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FORMULIZE0 for r/r c < 4: 



pi = 0.0818 + 0.0642/3 2 \J^-J + 0.0469/3 ^ 
+0.00626 (j-^j + 0.00393/3 2 f ^ 
-0.0317/3 f ^ ; 

p 2 = 0.671 + 0.0389 f^-\ + 0.04 1 8/3 2 f ^- 



-o.j in ; - 0.07SO ( I-j - 0.318,8 f^- 



-0.00495 — 
ft = 1.69 + 0.523 + 0.0130 f^- 

-0.417,0 - 0.717 - 0.140 f^- 



(13) 



(14) 



(15) 




200 300 400 
Index of profiles 

Figure 5. \ 2 normalized by the number of bins for the profiics 
in our test run. 



3.3 Parameter Space Test 

In obtaining the results in the above, we have made several 
simplifications. For example, the photons are collected from 
rings of finite width; particular values of £, T, k c and n s were 
adopted. We argued that our result are nonetheless general, 
because to is the only important factor governing the modi- 
fication of the line profile, and varying a at fixed £ is equiv- 
alent to varying to; for a fixed ion profile, T only changes 
a, which we directly constrain; and we showed in H3.1l that 
for k c > fc c ,crit, our results are independent of the values of 
fc c and n s . Nonetheless, these assumptions are worth check- 
ing. To confirm their validity and further test our model, 
we compare our model against an ensemble of 735 simu- 
lated profiles. Parameters are randomly drawn from reason- 
able ranges: k c G [27r/r c , 20/r c ], n s G [-2,-1.2], C G [2,6], 
T G [2, 8] keV, /3 G [0.5, 1.0] and v rms G [0, 500] km/s. Each 
simulation produces 2 x 10 photons, which are collected over 
a set of rings at r < 7r c (the largest radius within which our 
model has been calibrated). The rings have random widths 
between 0.25r c and 4r c , a range roughly corresponding to 
the half-power-diameter (HPD) of Astro-H relat i ve to th e r c 
of nearby clu sters (see Table 2 of lShang fc Obi l)2012al ) for 
the HPD and lChen et~afl l |2007f ) for the r c ). 

Fig. [5] shows the ratio of \ 2 to the number of bins rib 
for all the profiles in the test run. This ratio, similar to the 
reduced y 2 , is an indication of the goodness of fit. All of 
our profiles have X%/ n b ~ 1, except for a few outliers which 
we shall discuss later. The constraints on turbulence come 
primarily from the the height and variance of the line. It is 
therefore crucial that our model accurately reproduce these 
two quantities. In Fig. [6] we show the flux ratio between lines 
with and without RS. The x- and y- axes are the values from 
our model and the simulation, respectively, while the dotted 
line is where the model exactly agrees with the simulation. 
We see that resonant scattering can both decrease the flux, 
r < 1, by scattering photons from the cluster core out of the 



http:/ /formulizc. nutonian.com/ 




'flux, mod 



Figure 6. Flux ratios from the simuiation (y axis) and our model 
(x axis) between iine profiles with and without RS. The straight 
dotted line is r fluXjSim = r flux ^ mod 



field of view, and increase the flux, r > 1, by scattering pho- 
tons from the core into the field of view. Overall, our model 
shows very good consistency with the simulation, except for 
a few outliers. Similar consistency is seen in Fig.[7]where we 
plot the variance ratio between lines with and without RS. 

As for the a few outliers in Fig.[5]and[6] careful examina- 
tion reveals that they are all related to emission at large radii 
(r ;>6r c ). This likely indicates that our calibration at large 
r is inaccurate. Two reasons might be responsible for this. 
First, due to the declining surface brightness, the photon 
counts at large radii in our simulation are not large enough 
to give us good statistics. Second, the large photon collect- 
ing ring picks up residual large scale motions in the velocity 
field. Note that the velocity power spectrum did not go to 
as k — > 0, but stays as a constant. The parameter values 
in Table IA1I also show large scatter at large r. Our results 
for r ,>6r c therefore cannot be trusted; however, this region 
contributes negligibly to an emission-weighted observation, 
and this shortcoming is of little importance. 



© 0000 RAS, MNRAS 000, 000-000 



8 C. Shang and S. P. Oh 





1.4 r 




1.35 - 




1.3 - 




1.25 - 


E 


1.2 - 


<n 


1.15 - 


ca 




1.1 - 




1.05 - 




1 - 




0.95 - 




0.9 

0.9 0.95 1 1.05 1.1 1.15 1.2 1.25 1.3 1.35 1.4 

'Var.mod 

Figure 7. Variance ratios between line profiles with and without 
RS. As in Fig. [6] x- and y- axes are the results from our model and 
the simulation, respectively, while the dotted line is rg ux s j m = 

^flux,mod- 
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Figure 9. Electron (dashed green curves) and Fe XXV ion den- 
sities (red squares) as a function of radius for four nearby bright 
clusters. The dotted blue curves are the best-fit beta model to 
the ion profile outside the cores. 
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Figure 8. Comparison between our model and simulations, when 
the core radii for electrons and ions are the same, but the outer 
profiles differ. The ion density is assumed to follow a model with 
jS = 0.8, while the electron distribution has = 0.6. The dotted 
blue curves are model outputs when the ion profile is used as our 
model input while dashed green curves use the electron profile. 
The dot-dot-dashed curves are simulated PDFs without RS, while 
the data points are those with RS. The left and right panels are 
for the inner 0.5r c and an annulus from r c to 2r c , respectively. 
Using ion profiles provides accurate results which agree with the 
full simulations. 



3.4 Application to Clusters 

In the above, we performed calibrations for the case when 
the ions and electrons have identical ft profiles which differ 
only in normalization (which would require isothermal clus- 
ters of constant metallicity). However, distributions in real 
clusters deviate from such a idealized model, and will affect 
where the emissivity profile (oc n e rii), which dictates where 
photons originate, and the optical depth profile (oc m), 
which dictates where they scatter. We now discuss how these 
deviations affect our results, and how they can be corrected 
for. 

Most non cool-core clusters (NCCCs) can be rea- 
sonably be described by our model. This may seem 
somewhat surprising given that both of their tempera- 
ture and metallicity profiles tend to decline toward large 
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Figure 10. Comparison of line profiles between our model and 
simulations for four nearby CCCs. The line and point types are 
same as in Fig. [8] The PDFs are from the central 30 kpc. 



radii (e.g., |Pe Grandi fe Molendil 1200 ll ; iBaldi et al.l 120071 : 
iLeccardi. Rossetti fe Molendil 120101 )1 However, the abun- 
dance of Fe XXV, the ion responsible for the He-like iron 
line, declines with temperature at T >~ 3 keV (which is 
usually true for bright clusters). So the decline of temper- 
ature and decline of metallicity have opposite effects which 
roughly cancel, making the ion abundance roughly constant 
in the core. 

To conservatively estimate the uncertainty introduced 
by differences in the electron and ion profiles outside the 
core, we perform the following illustrative calculation. We 
assume the ion density and electron density are described 
by two P mo dels with the same core radius but different /3. 
According to lBaldi et al.l l|2007l ). both the temperature and 
metallicity decline with a power law slope of ~ 0.3 outside 
the core. Since n oc r~ 313 for r 2> r c , we set 3(/3i — f3 c = 0.6 
to account for this, or /3i — /3 C = 0.2. We then ran Monte 
Carlo radiative transfer calculations with these electron and 
ion profiles, and compared them to results from our fitting 
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Figure 11. Same as Fig, 1101 except that the PDFs here are ex- 
tracted from an annulus between 30 — 60 kpc. The thick dotted 
blue curves are when the "point source" correction is ap- 
plied (see text for details). 



formulae. The results are shown in Fig. [8] The left and right 
panels show the PDFs from the inner 0.5r c and an annulus 
from r c to 2r c , respectively. The simulated PDFs (the data 
points) are compared to our model with two sets of input 
parameters. The dashed green (dotted blue) curves are when 
the P model parameters for the electron (ion) distribution 
are used. The fit is good when the beta model parameters 
for the ion distribution, rather than the electron distribu- 
tion, is used. This is reasonable, since the region outside the 
core is more important for scattering rather than producing 
photons, so here the optical depth profile is more important 
than the emissivity profile. 

However, the brightest nearby clusters which Astro-H 
will first turn to are cool-core clusters (CCCs). In CCCs, 
a P profile can be a relatively poor fit to the ion distribu- 
tion in the center. In particular, observed temperature and 
metallicity profiles often imply ion profiles which are non- 
monotonic and can decline, rather than remain constant, 
toward the center. This deviation will affect the number of 
scatterings taking place in the core and the number of pho- 
tons emitted by the core but scattered at outskirts, affecting 
the profiles seen in the core and outskirts respectively. To 
estimate its importance, we perform calculations for four 
nearby bright CCCs - 2A0335, A0496, A3526, and Perseus 
- with measured temperature and density distributions. The 
ele ctron density distr ibutions are taken from the /3 model fit 
in IChen et al.l \2QQ% , whil e the temperature an d metalli- 
cally distributions are from lSnowden et al.l (|2008l ). Three of 
the four clusters, except A0496, show dips in their metal- 
licity distribution in the innermost bins, which might have 
been caused by RS. We thus ignore these bins, and instead 
assume the metallicity distribution is flat within the second 
innermost bin. We then compute their Fe XXV ion density 
profiles, based on these distributio ns and appropriate rate 
coeffi cients for coronal equilibrium IjBryans. Landi fe Savinl 
2009). They are shown in Fig. [9] with squares. Compared to 
the electron density distributions (the dashed curves), the 
ion distribution is steeper outside the cores. Notably, the ion 
distribution in A3526 has a steep drop in the center, due to 



the temperature drop. Overall, the complex ion distributions 
cannot be fully described by a simple and general form. 

Nonetheless, to test the applicability of our previous cal- 
culations, we perform /3 model fits to the profiles at r > 15 
kpc, shown with dotted blue curves. Monte-Carlo radiative 
transfer simulations were then performed for all these clus- 
ters, with line emissivities and oscillator strengths taken 
from the Astrophysical Plasma Emission C ode (APEC) and 
ATOMDE03 database l|Foster et alj|2012l ). The rms turbu- 
lent velocity is set to Oturb = 100 km/s. We compute the the 
PDFs from the central 30 kpc, corresponding to 0.5 -1 r c and 
spanning angles of 0.73 - 2.3 arcmin for these clusters. The 
spanning angles are comparable to the HPD ~ 1.3 arcmin of 
Astro-H. Fig [10] shows the comparisons between our model 
and simulations, where the dashed and dotted curves are 
when the electron and ion distributions are used as inputs. 
Our model is in good agreement with simulations when the 
P model of the ion distribution is used as the input. The 
flux and variance ratios agree within 8% and 4%, respec- 
tively. As expected, the agreement becomes worse when the 
electron distribution is used as the input. In Fig. 111! we com- 
pare PDFs accumulated between 30-60 kpc. The discrepancy 
between the model and simulation is somewhat larger, but 
still perfectly acceptable. The prominent exception is A3526, 
where the ion and electron profiles differ significantly in the 
center. Here, the PDFs are significantly overestimated by 
our model in these outer regions if the P model fit of the ion 
distribution is used as the input. These results can be under- 
stood as follows. RS affects the cluster center primarily as a 
source of effective absorption, by scattering photons out of 
the line of sight. All that matters in this case is the central 
line optical depth (as we shall see in i]3.6|l . which is set by 
the ion profile. However, further out, RS is also important as 
a source of re-emission, by scattering photons emitted from 
the cluster core back into the line of sight. While adopting 
the cuspy ion profile here correctly estimates r(r), it incor- 
rectly overestimates emission from the core oc n;n e , thus 
predicting an overabundance of scattered photons. We find 
that a simple scheme could be used to largely correct this 
bias. With the measured temperature and metallicity pro- 
file, we first compute the overestimated emission from the 
core. Then, assuming this emission is from a point source 
in the cluster center, we subtract their contributions from 
the line profile. The corrected line profiles, shown in Fig. [10] 
with thick dotted blue curves, follow the simulated profiles 
more closely. In particular, the agreement is now excellent 
for A3526. In order to apply this correction, a new calibra- 
tion has been performed for a point source in the cluster 
center. For this calibration, we use our parametric form to 
fit line profiles enclosed in a radius r, instead of those from 
a differential area at r. The latter could not be well- fitted by 
our model, since the emission come from scattered photons. 
To apply the model, the original emission from the point 
source has to be included. The profile from an infinitesi- 
mal area should be obtained from the calibration by taking 
a derivative dP/(2nr)dr. This "point-source" calibration is 
also publicly available on the internet. 

Based on these results, we conclude that our fits are 
broadly applicable to both CCCs and NCCCs. Using the 



http:/ /www. atomdb.org/indcx.php 



© 0000 RAS, MNRAS 000, 000-000 



10 C. Shang and S. P. Oh 



ion ft profiles generally provides more accurate results than 
using the electron /3 profiles. Where they are observed to 
differ significantly (as in A3526), a "point source" correction 
should be performed to correct line profiles originating from 
outside the core. 



3.5 Error Estimates and Analysis of Full Line 
Complex 

In this section, we discuss statistical and systematic errors 
in the moments of the line profile introduced by RS. By ap- 
plying our model to the full 6.7 keV Fe line complex, which 
harbors both optically thick and thin lines, we show that 
RS has a fairly mild effect on statistical uncertainty, but a 
large effect on systematic errors. Note that some fraction 
of the latter carries over if RS is imperfectly modeled. In 
the Gaussian model, the constraint on turbulence are solely 
driven by the variance of the line. RS modifies the line pro- 
file, changing both the height and its shape. So there are 
three sources of constraints: the line dispersion, the total flux 
in the line, as well as deviations from Gaussianity. Although 
the last potentially c ontains extremely valuable information 
l|Shang fe Ohll2012bl ). the effects of RS can be more compli- 
cated there, and we defer a detailed investigation to future 
work. We therefore focus on the first two constraints, which 
correspond to the zeroth and second moments of the line 
profile. 

It is easy to see why resonant scattering causes a mild 
increase in the statistical uncertainty of turbulent line broad- 
ening: it increases the dispersion of the line, and decreases 
the flux (or the number of photons) in a line. These effects 
are shows in Fig. 1121 where different line types correspond 
to different radii and width of photon collecting rings. Flux 
is mostly absorbed in the inner regions rfl ux < 1, and re- 
emitted in the outer regions rfl ux > 1. Note that the differ- 
ence in dispersion between profiles with and without reso- 
nant scattering, 8a — a-R,s — a OI i ~ const, independent of 
cr or i. This allows us to write: 



AlJtur 



AfTturb.ori 



1 + 



da 



(16) 



which is greater than unity, since rfl ux < 1 and 5a > 0. 
However, as we shall see, we expect errors to be dominated 
by systematic errors (due to imperfect modeling of resonant 
scattering, which has a large effect), rather than statistical 
errors (which are small, due to the large number of photons 
n ph ~ 10 4 . 

We now perform a fully realistic spectral line fit of the 
entire 6.7 keV Fe line complex. Besides allowing a fully ac- 
counting of the error budget, it highlights the systematic 
biases introduced by RS if it is ignored. It also serves as a 
final check of our model. We first produce mock data of line 
profiles from the inner 30 kpc of the Perseus cluster, using 
the temperature, density and metallicity profiles described 
in §E1 The ID RMS velocity is set to a tulh = 100 km/0. 
We renormalize the profiles such that the photon counts in 
the Helium-like iron line is 10 4 , roughl y the expected count s 
of Astro-H from local bright clusters jShang fe Ohll2012ah . 



11 This number is illustrative; there are hints that turbulence in 
Perseus could be substantially stronger iChurazov et alj|2004l . 
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Figure 12. Dispersion difference 8a = .(Jus ~ fori and flux ratio 
as functions of <T OI i, the line width without RS. These lines all as- 
sume C = 5 an d P = 0-6- Photon are accumulated over rings with 
{inner radius, width}={0, 0.5r c } (solid red line), {0, r c } (dashed 
green line), {0, 1.5r c } (dotted blue line), {r c ,0.5r c } (dot-dashed 
purple line) and {3r c ,0.5r c } (dot-dot-dashed brown line). 



The counts in other lines and in the continuum are assigned 
according to the APEC model l|Foster et al.ll2012T ). assum- 
ing a metallicity of 0.5 Zq. Convolving with the telescope 
response (assuming FWHM = 5 eV), we produce mock spec- 
tra both with and without RS. Most lines in the complex 
usually have very small optical depth, except the Helium-like 
iron line. For simplicity, we here assume only the Helium-like 
iron line is affected by the RS. We then use the Gaussian 
and the new model to fit the spectra, adopting \ 2 statistics 
(equation [TTJ and the Markov chain Monte Carlo al gorithm 
implemented in CosmoMC (|Lewis fe Bridle! 120021 ') . There 
are four free parameters: the original photon counts in the 
Helium-like iron line n p h, the mean [i, the original width a 
and the metallicity Z. We studied 3 cases: the new model fit- 
ting the RS spectrum, the Gaussian model fitting the Gaus- 
sian spectrum (unmodified by RS) and the Gaussian model 
fitting the RS spectrum. The results are given in Table Q] 
Fig. [13] shows the marginalized posterior distributions. They 
show our model correctly reproduces these quantities. The 
constraints are a little worse than those in the case 1, as 
expected. In case 2 where the effects of RS are ignored, the 
parameters are significantly biased. The line width is over- 
estimated due to the additional broadening caused by RS. 
The metallicity is on the other hand underestimated due to 
the reduced counts in the line. An improved model like ours 
is therefore imperative to correctly extract the underlying 
parameters. 
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Figure 13. ID marginalized posterior distributions for /x, a and 
Z. Solid red, dashed green and dotted blue curves correspond to 
case 0, 1 and 2, respectively. The vertical black lines indicates the 
location of the input values. Ignoring resonant scattering results in 
significant biases in the recovered dispersions a and metallicities 



Table 1. Best- fit values and 68% uncertainties from fitting the 
mock spectra. Case is when our model is applied to a RS spec- 
trum, case 1 is when the Gaussian model is applied to a Gaussian 
spectrum (i.e., without RS), and case 2 is when the Gaussian 
model is applied to a RS spectrum. 



n ph (xfO 4 ) ,u(eV) 



a (eV) 



Z (Z e ) 



Input 
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2.82 



0.3 



line photons in moderately optically thick media, where the 
natural line wings are unimportant and only the Doppler 
core has significant optical depth. In this case, for line center 
optical depths To, for an initial line profile J(x) — je~ x /i/tt 
(i.e., x = AE/V2a in our notation), the final escaping lin e 
profile for isotropic scattering with no recoil is (|Fieldlll959l ): 



-J(x) _ 1 — cxp(— T e x ) 
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where 



9n(U) 
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J )du 



(17) 
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(log W -l)l/2 JQ 

We have found that this integral can be evaluated analyti- 
cally: 
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where y = (logu -1 ) 1 / 2 . 

This profile was derived by considering the time devel- 
opment of the frequency distribution o f a photon em itted in 
a nebula of infinite optical depth, but iFieldl (Il959l ) showed 
that it can also be applied with reasonable approximation to 
a nebula with finite optical depth to. Focusing on the first 
term, we see that there is a critical frequency: 

t„ = 1 => x c = (logr ) 1/2 (20) 

at which the wings of the profile become optically thin, and 
photons can escape freely. Indeed, for x > x c , the derivative 
of the profile is approximately — 2xe~ x l\pn, its initial value, 
while for x < x c , the derivative is exp(— Toe~ x ), and hence 
quickly approaches zero. Thus, the profile is flat for x < x c , 
and relatively unchanged from the initial profile for x > x c 
— exactly what we have seen in our Monte-Carlo simula- 
tions. The escape probability per scattering is erfc(Vhi to), 
and the mean number of scatterings is (N) = l/erfc(\/hiTo). 

Equation [17] was obtained by explicit solution of the 
integro-differential equation of radiative transfer. Nonethe- 
less, we can intuitively understand this line shape from the 

= form of the redistribution function r(x,x'), which gives the 
probability that radiation of frequency x will be s cattered 

" into frequency x ' . For isotropic scattering, it is l|Zanstral 
ll949l : lFieldlll959l ): 



, f°° _„2 

r(x, x') = I e " dv; \ 



x\ = max(|a;|, |a;'|). (21) 



3.6 Analytic Model for the Effects of RS 

In this section, we discuss an approximate analytic model 
for the modifications in the line profile due to resonant 
scattering. This drives both physical understanding, and al- 
lows very rapid estimates, particularly of important physical 
quantities such as modifications to the line height, variance, 
and flux. 

The previously calculated line profiles in the central re- 
gions have a characteristic shape: flat in the central regions, 
matching onto Doppler wings with relatively little change 
from the unscattered profile. This characteristic shape is a 
generic feature arising from radiative transfer of resonance 



-where v is in units of the ID velocity dispersion a v ; note 
that r(x,x') is symmetric. Physically, we can understand 
this expression from the fact that the speed of an atom must 
be larger than both x and x' to scatter x into x'; as long 
as this is true, then for isotropic scattering the scattering 
probability is independent of x and x' . Thus, for a photon 
of frequency x, the redistribution function is constant (inde- 
pendent of x'), until x' — x, at which point r(x,x') plum- 
mets rapidly. We can apply this to understand the form of 
the line profile. In the optically thick core x < x c , photons 
are trapped and scatter repeatedly. Due to the form of the 
redistribution function, the scattered photons have a flat 
profile between —x and x, with a finite chance to escape to 
x > x. One can think of the photons as being trapped be- 
tween reflecting boundaries at —x and x, which have some 
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finite probability of transmission. Thus, the Doppler profile 
in the core is erased to form a flat profile. Photons gener- 
ally escape only once they hit a rare fast moving atom, and 
are re-emitted with x > x c . Photons initially emitted in the 
optically thin wings x > x c of course travel unimpeded to- 
wards us, and the line wings are unaltered from the initial 
profile. 

Besides redistribution in frequency, RS also causes spa- 
tial diffusion: photons migrate from the optically thick clus- 
ter core to the optically thin cluster outskirts, causing the 
surface brightness profile to be more heavily weighted to- 
ward the optically thin 'photosphere' than the emission pro- 
file (in the absence of any destruction processes, the total 
line flux from the cluster is, of course, conserved; the to- 
tal area under the curve in equation (|17p is unity). This 
can be clearly seen in Fig 1121 where the flux ratio in the 
line with and without resonance scattering is rfl ux < 1 in 
the inner regions and rfl ux > 1 in the outer regions. Also, 
from Table I All the parameter g\ , which parametrizes the 
strength of re-emission, increases with radius. We therefore 
have to be careful when comparing the analytic formula in 
equation (|17[) to simulation results. Let us focus on the fate 
of photons which are generated in the inner cluster core; 
the fraction of photons which originate in the cluster out- 
skirts is small. When considering the emission from radial 
line profiles P(r), in the inner regions the photons which 
are scattered out of the line of sight (to be re-emitted in the 
outer regions) are effectively 'absorbed', since they do not 
contribute to emission from that radius. Thus, only the first 
term in equation \17\j should be used. However, when con- 
sidering the observed line profile emerging from the cluster 
as a whole (which involves an integral over all radii), both 
terms in equation (|17[l should be used, since in this case the 
photon flux is conserved. 

This analytic formula requires the optical depth at line 
center. Note that this is a spatially varying quantity, since 
a photon from the near side of a cluster scatters less than a 
photon from the far side of the cluster. For a shell at radius 
r, the correct optical depth to use is: 

reff(r) = -log (^JdnJ^ dl exp[-n(n,r)al]\ (22) 

= -log Q J dcos{9) J dl exp[-n(\/' 2 + r 2 - 2rlcos{e))ol] S j 

For shells close to the cluster center, the correction is small, 
and a simple integral from r to infinity suffices. 

Fig. [1] compares the analytic profile to simulation re- 
sults; they show good agreement. An analytic profile allows 
simple statistics such as modifications to line height, vari- 
ance and flux due to resonant scattering to be quickly esti- 
mated; we find that these reproduce the trends seen in Fig 
1121 and are accurate to ~ 20 — 30%. 

We stress that the analytic profile is only approximate; 
quantitative conclusions should always been drawn from ra- 
diative transfer simulations. The analytic profile is accurate 
near the cluster center and for the spectrum emerging from 
the cluster as a whole, but modeling the radial trends seen in 
our simulations (which reflect photons from the cluster core 
scattered into the line of sight) is beyond the scope of this 
paper. The analytic approach also makes certain restrictive 
assumptions, though these turn out not to matter signifi- 
cantly. Equation (|17|) is derived under the assumption that 



ro 3> 1, though in practice it should be reasonably accurate 
for moderate to > 1. Note also that equation (|17|l also as- 
sumes isotropic scattering, whereas for instance our fiducial 
6.7 keV He-like Fe is pure dipole scattering. Dipole scat- 
tering has a different redistribution function r(x,x') from 
isotropic scattering; there is a correlation between the in- 
going and outgoing photon frequencies, since the scattering 
tends to be in the same direction as absorption. In prac- 
tice, we find from our Monte-Carlo simulations that calcu- 
lations which assume isotropic and dipole scattering do not 
differ significantly. Overall, the analytic approach is useful 
for physical insight and quick 'rule of thumb' estimates. 



4 SUMMARY 

Upcoming telescopes will enable us, for the first time, to 
measure ICM turbulence through Doppler broadening of 
emission lines. However, the brightest lines with the highest 
signal to noise also tend to be mildly optically thick. Res- 
onance scattering of emission lines needs to be taken into 
account, as it alters the line flux and more importantly, in- 
creases the line widths used to measure turbulence. Naive 
measurements which ignore this would give biased results. 
To date, characterizing and eliminating this bias has not 
been quantitatively addressed. In this paper, we proposed a 
new model for the emission line profile, accounting for the 
effects of resonant scattering. This model assumes that the 
ions and electrons both obey /3 profiles, and works best in 
the limit where these profiles are fairly similar in shape in 
the core. 

Motivated by cosmological simulations of galaxy clus- 
ters, we assume that the velocity field in galaxy clusters is 
adequately characterized by a power spectrum alone, and 
show that as long as the outer scale is smaller than the 
core size, for a given rms velocity, our results are not sensi- 
tive to the assumed power spectrum. We run Monte Carlo 
radiative transfer simulations and show that our model pro- 
vides an excellent fit with only 3 parameters, allowing the 
effects of resonant scattering and turbulent broadening to 
be disentangled. We provide fitting formulae and tables for 
model parameters which will allow rapid fits to observed 
line profiles. As a test, we run MCMC simulations of the 
entire Fe line complex for the Perseus cluster, and show our 
model provides unbiased estimates of turbulent broadening 
with relatively small error bars, where neglecting resonant 
scattering clearly results in biased estimates of turbulent 
widths and metallicities. Finally, we show that the physical 
origin of the calculated line profiles — which are typically 
flat in a central core and have relatively unmodified wings 
- can be understood analytically. Our model assumes that 
the electron profile (which affects the emissivity profile) and 
the ion profile (which affects both the emissivity and opti- 
cal depth profile) have relatively similar /3 profile shapes in 
the core; where they differ, the ion profile should be used. 
When they differ significantly, a "point source" correction 
should be used to correct emission originating from outside 
the core. In this way, CCCs are well fitted by our model as 
well as NCCCs. 

Several extensions and improvements to this work are 
possible. Our discussions have focussed on the He-like 6.7 
keV Fe line; detailed considerations of other lines would 
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be wor thwhile. Secondly, as we discussed in IShang & Ohl 
l|2012al ) , distinctive large scale modes of gas motions arise in 
clusters. These modes could potentially be separated and in- 
dividually constrained by mixture modeling. Such mode sep- 
aration could be very useful separating bulk motions (e.g., 
gas sloshing in clusters with cold fronts) from turbulence, 
and discerning volume filling fraction of turbulence. The ef- 
fect of resonant scattering has yet to be taken into account 
in mixture models, which exploit non-Gaussian features in 
the line profile. 
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APPENDIX A: BEST-FIT MODEL 
PARAMETERS 



Table Al. Best-fit model parameters for different and r. 



r/r c 




PI 






P2 






h 




= 0.5 


0. 


05 


0. 


.098 ± 





.027 


0.615 ±0.028 


1 


460 


±0 


.224 





.55 


0. 


092 ± 


0. 


.008 


0.489 ± 0.009 


1 


.359 


±0 


.074 


1 


05 


0. 


122 ± 





.007 


0.396 ± 0.007 


1 


.089 


±0 


.051 


1. 


.55 


0. 


124 ± 





008 


0.301 ±0.007 


1 


100 


±0 


.049 


2 


05 


0. 


141 ± 





008 


0.237 ± 0.008 


1 


.012 


±0 


.047 


2 


.55 


0. 


153 ± 





009 


0.185 ±0.008 





.955 


±0 


.047 


3 


.05 


0. 


154 ± 





.009 


0.140 ±0.009 





.942 


±0 


.049 


3 


.55 


0. 


.153 ± 





010 


0.105 ±0.009 


0. 


.921 


±0 


.053 


4. 


05 


0. 


165 ± 





010 


0.082 ± 0.009 


0. 


.862 


±0 


.054 


4. 


.55 


0. 


169 ± 


0. 


Oil 


0.067 ±0.010 





.844 


±0 


.056 


5 


05 


0. 


164 ± 





.012 


0.048 ±0.011 





933 


±0 


.060 


5 


.55 


0. 


167 ± 





.012 


0.035 ±0.011 





.912 


±0 


.061 


6 


.05 


0. 


140 ± 





013 


-0.005 ±0.012 





.907 


±0 


.077 


6 


.55 


0. 


153 ± 





014 


0.002 ±0.013 


0. 


.922 


±0 


.076 


7. 


.05 


0. 


164 ± 





.015 


0.006 ±0.013 


0. 


.938 


±0 


.074 


= 0.6 





.05 





074 ± 





019 


0.550 ±0.021 


1 


.541 


±0 


.213 





.55 


0. 


.095 ± 





.007 


0.464 ± 0.007 


1 


.274 


±0 


.056 


1 


.05 


0. 


128 ± 





.006 


0.352 ± 0.006 


1 


.103 


±0 


.040 


1 


.55 


0. 


147 ± 





.007 


0.256 ± 0.007 


1 


.002 


±0 


.039 


2. 


05 


0. 


161 ± 





008 


0.178 ±0.007 


0. 


.946 


±0 


.041 


2 


.55 


0. 


.172 ± 





009 


0.123 ±0.008 


0. 


861 


±0 


.044 


3 


.05 


0. 


175 ± 





.010 


0.082 ±0.009 





.878 


±0 


.047 


3 


55 


0. 


.174 ± 





011 


0.051 ±0.010 





.850 


±0 


.053 


4. 


05 


0. 


165 ± 





012 


0.012 ±0.011 


0. 


.856 


±0 


.062 


4. 


.55 


0. 


169 ± 


0. 


.013 


-0.006 ±0.012 





.863 


±0 


.065 


5 


.05 


0. 


165 ± 





.014 


-0.030 ±0.013 





.874 


±0 


.071 


5 


.55 


0. 


.154 ± 





015 


-0.047 ±0.014 


0. 


874 


±0 


.082 


6 


.05 


0. 


.173 ± 





016 


-0.044 ±0.015 


0. 


.820 


±0 


.081 


6 


.55 


0. 


.174 ± 





018 


-0.040 ±0.016 


0. 


.931 


±0 


.085 


7. 


.05 


0. 


182 ± 





.019 


-0.040 ±0.017 





.860 


±0 


.090 


= 0.7 


0. 


05 





088 ± 





015 


0.545 ±0.016 


1 


411 


±0 


145 





.55 


0. 


100 ± 





.005 


0.439 ±0.006 


1 


.256 


±0 


.045 


1 


05 


0. 


140 ± 





.006 


0.308 ± 0.005 


1 


.028 


±0 


.033 


1 


55 


0. 


166 ± 


0. 


.006 


0.197 ±0.006 





.904 


±0 


.034 


2 


05 


0. 


175 ± 





008 


0.114 ±0.007 


0. 


.886 


±0 


.038 


2. 


.55 


0. 


.179 ± 





009 


0.054 ±0.008 


0. 


.906 


±0 


.044 


3 


.05 


0. 


195 ± 





Oil 


0.025 ±0.009 


0. 


.855 


±0 


.048 


3 


55 





203 ± 





012 


-0.003 ±0.011 


0. 


.826 


±0 


.053 


4. 


05 


0. 


210 ± 





.014 


-0.020 ±0.012 


0. 


.852 


±0 


.059 


4. 


.55 


0. 


210 ± 


0. 


.016 


-0.048 ±0.014 





.791 


±0 


.069 


5 


.05 


0. 


206 ± 





.018 


-0.060 ±0.015 





.833 


±0 


.075 


5 


,55 


0. 


.204 ± 





020 


-0.082 ±0.017 


0. 


851 


±0 


.085 


6 


.05 





208 ± 





022 


-0.083 ±0.019 





783 


±0 


.093 


6 


.55 


0. 


.205 ± 





.025 


-0.091 ±0.022 





.880 


±0 


103 


7. 


.05 


0. 


240 ± 


0. 


.028 


-0.085 ±0.024 





.742 


±0 


.106 
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Table A2. Continued from Table [ATI 



r/r c 




Pi 




P2 






h 




P = 0.8 





.05 


0. 


.105 ±0. 


.013 


0.540 ±0.013 


1. 


.277 


±0. 


101 





.55 


0. 


113 ±0. 


.005 


0.416 ± 0.005 


1. 


.218 


±0. 


035 


1. 


.05 


0. 


156 ± 0. 


.005 


0.267 ±0.005 


0. 


.995 


±0. 


028 


1 


.55 


0. 


173 ± 0. 


.006 


0.139 ± 0.006 





.889 


±0. 


032 


2. 


.05 


0. 


196 ± 0. 


.008 


0.058 ± 0.007 


0. 


.808 


±0. 


.037 


2 


.55 


0. 


203 ± 0. 


.010 


-0.005 ± 0.009 


0. 


.816 


±0. 


.046 


3 


.05 


0. 


236 ± 0. 


.013 


-0.022 ±0.011 





.807 


±0. 


049 


3 


.55 


0. 


240 ± 0. 


.015 


-0.062 ± 0.013 





.772 


±0. 


059 


4. 


.05 


0. 


235 ± 0. 


.018 


-0.097 ±0.015 





.737 


±0. 


070 


4. 


.55 


0. 


271 ± 0. 


.021 


-0.086 ±0.017 


0. 


.762 


±0. 


071 


5. 


.05 


0. 


258 ± 0. 


.025 


-0.122 ±0.020 





.821 


±0. 


085 


5. 


.55 


0. 


340 ± 0. 


.030 


-0.065 ± 0.022 





.827 


±0. 


.078 


6. 


.05 


0. 


309 ± 0. 


.035 


-0.111 ±0.027 


0. 


.818 


±0. 


099 


6. 


.55 


0. 


323 ± 0. 


.040 


-0.115 ±0.031 


0. 


.856 


±0. 


104 


7. 


.05 


0. 


391 ±0. 


.048 


-0.098 ± 0.034 


0. 


.767 


±0. 


109 


P = 0.9 





.05 


0. 


.099 ± 0. 


Oil 


0.521 ± 0.012 


1 


.325 


±0. 


094 





.55 


0. 


120 ± 0. 


.004 


0.389 ± 0.004 


1 


.164 


±0. 


.030 


1. 


.05 


0. 


162 ± 0. 


.005 


0.221 ± 0.005 





.972 


±0. 


026 


1. 


.55 


0. 


194 ±0. 


.007 


0.099 ± 0.006 


0. 


.862 


±0. 


030 


2. 


.05 


0. 


217 ±0. 


.009 


0.011 ± 0.008 





.812 


±0. 


038 


2. 


.55 


0. 


254 ± 0. 


.012 


-0.040 ± 0.010 


0. 


.715 


±0. 


.045 


3 


.05 


0. 


286 ± 0. 


.016 


-0.058 ±0.013 


0. 


.821 


±0. 


.050 


3 


.55 


0. 


324 ± 0. 


.020 


-0.079 ± 0.016 





.757 


±0. 


058 


4. 


.05 


0. 


330 ± 0. 


.025 


-0.121 ±0.019 





.713 


±0. 


.072 


4. 


.55 


0. 


346 ± 0. 


.032 


-0.141 ±0.024 


0. 


.795 


±0. 


084 


5. 


.05 


0. 


367 ±0. 


.039 


-0.137 ±0.028 


0. 
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